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Symbols 


d 

mean temperature' difference 4 , K, (see eq. (1)) 

m 

sample size 

V 

fi 

pressure?. hPa (or 100 Pa) 

sam])le standard deviation, K. (sec? eq. (2)) 

T 

temperature, K 

t 

Student’s /-test, statistic* (sec? eq. (3)) 

X 

station temperature value' from single' sounding. K 

Z 

geopotential height, km 

z 

geometric altitude, km 

Subscripts: 

J 

pair number (see eq. (1)) 

m 

sample size (see eq. (3)) 

1 

satellite (see eq. (1)) 

2 

in situ (see eq. ( 1 )) 

Abbreviations: 


FSU 

former Soviet Union 

IFS 

inflatable falling sphere 

LAIPAT 

LIMS inverted profile arc hive tape? 

LAMAT 

LIMS map arc -hive tape? 

LIMS 

Limb Infrared Monitor of the Stratosphere’ 

MI) 

me'mi difference 

NMC 

National Meteorological Center 

NO A A 

National Oceanic and Atmospheric Administration 

RAOB 

radiosonde' observation 

ROCOB 

rocket sonde observation 

SAO 

semiannual oscillation 

SD 

standard deviation 

ssu 

stratospheric sounding unit 

TIROS 

Television and Infrared Operational Satellite' 

TOYS 

TIROS Operational Vertical Sounder 

UTC 

universal time coordinate 

VTPR 

Vertical 'temperature* Profile' Radiometer 

WMO 

World Meteorological Organization 
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Summary 

The Limb Infrared Monitor of the Stratosphere 
(LIMS) experiment made observations from the Nim- 
bus 7 satellite in 1978 and 1979. Temperature-versus- 
pressure, T(p), profiles were derived from its limb 
radiance measurements, and those profiles were used 
to register the radiances from other channels and to 
retrieve species concentration profiles from those ra- 
diances. Therefore, biases in the T(p) results must 
be known in order to estimate the accuracy of those 
species profiles. LIMS temperatures have been vali- 
dated in the past with colocated radiosonde and rock- 
etsonde measurements. The present report describes 
time series comparisons between satellite and rocket- 
sonde T(p) values at station locations. This approach 
to validation retains nearly all the rocketsonde pro- 
files, increasing sample size significantly (to 665). As 
a result, one can know better whether there is a bias 
that varies as a function of pressure altitude, latitude, 
or season. 

The results indicate no clearly significant bias for 
LIMS versus Datasonde from 10 hPa to 1 hPa at low 
and mid latitudes. There is a positive LIMS bias of 2 
to 3 K in the upper stratosphere at high latitudes 
for the Northern Hemisphere in both winter and 
spring. LIMS is progressively colder than Datasonde 
from 0.4 hPa (about -3 K) to 0.1 hPa (about —9 K) 
at all latitudes. A similar comparison between LIMS 
and the more accurate falling sphere measurements 
reveals an equivalent mid-latitude LIMS bias at 

0. 4 hPa but a much smaller bias at 0.1 hPa (—4.6 K). 
Because the biases do not vary noticeably with sea- 
son, it is concluded that they are not a function of 
atmospheric state. This result confirms the robust- 
ness of the LIMS temperature retrieval technique. 

LIMS comparisons with the Soviet M -100 rocket- 
sonde show significant biases in both the strato- 
sphere and the mesosphere; the Datasonde is consid- 
ered more useful as a validation standard. National 
Meteorological Center (NMC)/Datasonde mean dif- 
ferences are very similar to those for LIMS/Datasonde 
at 10 and 5 hPa. However, the quality of the NMC 
comparisons is reduced at 2, 1, and 0.4 hPa, primar- 
ily due to a lack of nadir radiance data from those 
levels during the 1978-79 period. Standard devia- 
tions for the differences are generally larger at all 
levels for NMC data than for LIMS data, indicating 
that the LIMS analyses follow the true temperature 
variations better than the NMC analyses. 

1. Introduction 

Temperature-versus-pressure, T(p), profiles are a 
fundamental product of satellite midinfrared limb 


emission sounders of the middle atmosphere. Specif- 
ically, Gille and House (1971) and Bailey and Gille 
(1978) showed that one can retrieve T(p) by using 
observed radiances versus relative altitudes from two 
radiometer channels that view the atmospheric limb 
in the v<± (or 15 pm) region of the CO 2 spectrum. 
The “narrow” bandpass CO 2 channel, located near 
the center of that band, has an emissivity- versus- 
pressure profile that approaches a value of 1 near a 
“reference” altitude of 30 km. At that level the ef- 
fective radiating temperature is very close to the at- 
mospheric temperature, according to the blackbody 
function. That temperature is then used to calculate 
an effective emissivity in the more transparent (or 
“wide”) CO 2 channel. The “reference” pressure for 
the 30-km point is then determined from a curve of 
the “wide” channel emissivity versus pressure. The 
hydrostatic equation is then used to calculate the 
whole T(p) profile from the observed “wide” channel 
profile of radiance versus relative altitude. This pro- 
cess is iterated to achieve a final T(p) profile in the 
stratosphere. Absolute altitudes for each pressure 
level are determined later by a hydrostatic integra- 
tion of T(p) with the aid of an independent analysis 
of the height of the 50-hPa surface. 

The pressures associated with the T(p) results 
from the Nimbus 7 Limb Infrared Monitor of the 
Stratosphere (LIMS) experiment were used to reg- 
ister the measured radiance profiles from each of its 
channels (Gille and Russell 1984). Then, the T(p) 
values were applied to those radiances for the re- 
trieval of the concentration profiles of LIMS con- 
stituents O 3 , H 2 O, HNO 3 , and NO 2 . This means 
that the LIMS constituent concentrations can be af- 
fected by a bias in those temperatures. In particu- 
lar, whenever the tangent-layer signal becomes low, 
the retrieved concentrations are more sensitive to 
temperature bias. Error studies carried out on the 
LIMS species show that temperature bias error is 
the largest source of error in the species retrievals 
at most pressure levels. Because LIMS temperatures 
and species concentrations vary with pressure (or al- 
titude), latitude, and season, the effect of a temper- 
ature bias can change accordingly. 

Temperature observations from meteorological 
rocketsondes (ROCOB’s) and radiosondes (RAOB’s) 
are considered as correlative data for validating satel- 
lite measurements. In a comparison of colocated 
LIMS profiles with ROCOB’s and RAOB s, Gille 
et al. (1984a) showed that the mean differences were 
generally within ±2 K below the 1 -hPa level (altitude 
of about 48 km). The comparisons with ROCOB’s 
become less reliable at higher altitudes as sources 
of errors for ROCOB’s become more pronounced. 



Several rocket-sonde techniques have been employed 
in the past by different countries. To be able to de- 
rive maximum information from limited rocketsonde 
soundings, the compatibility of various systems was 
studied. Intercomparison campaigns were carried out 
at Kourou, French Guiana (Finger et al. 1975), and 
at Wallops Island, Virginia (Schmidlin et al. 1980). 
Gillc ct al. (1984b) also used LIMS temperatures as a 
transfer standard between ROCOB’s obtained with 
instruments from the US and the former Soviet Union 
(FSU)-- identified as “USSR rocketsondcs” in their 
paper and in the remainder of this report. 

Remsberg et al. (1984) also carried out temper- 
ature comparisons with correlative measurements as 
part of LIMS 0;j validation activities and found good 
agreement. Remsberg (1986) compared LIMS zonal- 
mean temperatures with a 4-year temperature cli- 
matology derived from Rayleigh backseat ter lidar 
measurements at 44° N during March, April, and 
May, when effects of zonal waves are weak. The 
agreement between those two data sets was better 
than 3.5 K between 37 and 64 km. 

Time series comparisons at single stations show 
differences in both the phase and the amplitude of 
temperature waves. For example, Miles et al. (1987) 
and Grose et al. (1988) carried out comparisons 
of LIMS temperatures with RAOB data at Inver- 
cargill, New Zealand, and Berlin, Germany, respec- 
tively, by employing the Fourier coefficient values 
on the LIMS map archive tape (LAMAT) product 
(Remsberg et al. 1990). In each case the phase of 
the LAMAT temperature time series agreed with the 
changes observed by the RAOB’s, but wfith some re- 
duction in temperature wave amplitude. Although 
the LIMS LAMAT product contains some spatial 
smoothing, one can determine a temperature value 
for the exact location of a correlative measurement 
station. Miles et al. (1987) reported LIMS LAMAT 
minus RAOB mean differences (MD) at 100 hPa that 
were loss t han 1.2 K with rrns differences of less than 
2.5 K. The agreement between RAOB and LIMS in- 
verted profile archive tape (LAIPAT) data at Inver- 
cargill was even better, with MD less than 0.3 K, and 
rms differences of 1.6 K. The LAIPAT comparisons 
w'erc limited to fewer samples, however. 

Recently, Remsberg et al. (1992) compared LIMS 
and National Meteorological Center (NMC) temper- 
atures with RAOB temperatures in the Arctic lower 
stratosphere. They used a set of 22 stations with 
nearly uniform longitudinal and latitudinal coverages 
from 60° N and 84° N. The LIMS and NMC temper- 
ature analyses were compared with RAOB data at 
pressure levels between 100 and 10 hPa. The LIMS 
temperatures showed very good agreement (MD less 


than ±0.3 K) with RAOB’s for pressure levels be- 
tween 70 to 30 liPa as averaged over a 7- month 
period. The corresponding NMC minus RAOB MD 
values were within ±0.3 K from 100 to 30 hPa. The 
satellite/ RAOB comparisons at 10 hPa were limited 
somewhat because of fewer radiosondes ascending to 
this level, especially during winter. Generally, the 
number of sondes that reach 10 hPa in winter is 
10 percent of the total (Mclnturff 1978). 

The time series comparison approach employed 
by Remsberg ct al. (1992) is very useful and is com- 
plementary to the colocated vertical profile compar- 
isons reported in Gille et al. (1984a). Any seasonal or 
shorter period component in one or the other dataset 
can be identified, in principle. Because of a decrease 
of R AOB data from altitudes at the 10-hPa level and 
above, ROCOB’s arc important for a validation of 
LIMS temperatures in the upper stratosphere and 
mesosphere. Therefore, in this report we extend the 
comparison process upwards in altitude, and present 
time series of both LIMS LAMAT and NMC tem- 
peratures versus US and USSR ROCOB’s between 
pressure altitudes of 10 and 0.1 hPa, as outlined in 
section 2. Parallel with the previous comparison ex- 
ercises, we then compare time series of station tem- 
peratures derived from the LIMS and NMC analy- 
ses. We consider factors that influence the accuracy 
of those data sets as well as the influence of their 
vertical resolutions. Furthermore, we report any lat- 
itudinal or seasonal trend between the data sets. 

ROCOB’s must be used with some caution. For 
example, Mclnturff (1978) reports that, on aver- 
age, about one-third of ROCOB’s were unusable for 
NMC’s weekly synoptic analyses in the 1970’s. While 
the acceptance rate for most stations was better than 
80 percent, only about 40 percent were “usable” from 
Primrose Lake, Canada, from Thule, Greenland, and 
from Thurnba, India. This relatively low' acceptance 
rate was not biased toward any particular sonde type. 

Recently, Schmidlin et al. (1991) showed that the 
Super Loki inflatable falling sphere (IFS) technique 
provides temperatures inferred from density determi- 
nations that are more accurate than the measured 
Super Loki Datasondc temperatures in the meso- 
sphere. The aerodynamic heating and radiative heat- 
ing/cooling effects, which can introduce large errors 
for Datasondc temperatures in the mesosphere that 
must be corrected, do not affect the inferred IFS 
temperatures. In the stratosphere, random-like, but 
small-scale, vertical structure is induced in the de- 
rived IFS temperature profile because of the effect of 
vertical winds on the lightweight sphere; the verti- 
cally smoothed temperature profile is more accurate. 
An absence of significant biases for the IFS technique 
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provides a means of standardization of Datasonde 
ROCOETs. By combining the results from the com- 
parison of the IFS and the US Datasonde temper- 
atures with our intercomparison studies, the accu- 
racy of the LIMS and NMC satellite temperatures is 
reassessed. 

Section 3 contains a brief description of the data 
sets. The comparison method and results are pre- 
sented for individual stations in section 4, and the 
findings for LIMS/ROCOB ’s arc reported by latitude 
zone in section 5. Conclusions about T(p) accuracy 
are discussed in section G. Two appendices present 
the data comparisons in detail, both graphically and 
in tabular form. 

2. Approach 

The LIMS/ROCOB comparison statistics in Gille 
ct al. (1984a) were obtained from individual LIMS 
profiles that met strict space and time colocation 
criteria with a rocket sounding. As a result, many 
rocket profiles were not included in their statistical 
sets. They did not report a monthly or even a sea- 
sonal statistical breakdown from their sample. Be- 
cause the present time scries comparison is based on 
the LIMS-mapped Fourier coefficient product, we are 
able to calculate a T(p) value at the exact station lo- 
cation from the coefficients. Although a mapping of 
the LIMS profiles leads to a smoothed temperature, 
we can include all the rocket profiles in our compar- 
isons and thereby increase sample size significantly. 
ROC OB comparisons with these smoothed satellite 
data lead to larger standard deviations (SD) for a 
set of paired observations, but the mean differences 
(MD) are not affected much (Miles et. al. 1987). A 
determination of changes of MD with pressure, lati- 
tude, or season is the primary goal of this study. 

The comparison results are presented in two ways. 
First, the LIMS/NMC/ROCOB time series plots are 
presented for each station at the standard NMC anal- 
ysis pressure levels of 10 hPa (approximate altitude 
31 km), 5 hPa (36 km), 2 hPa (43 km), 1 hPa 
(48 km), and 0.4 hPa (55 km). LIMS/ROCOB com- 
parisons are also reported at 0.1 hPa (65 km). Using 
the plots, one can visually evaluate whether the satel- 
lite data follow the observed temperature variations 
seen in a rocked, time series. In effect, the time se- 
ries plot represents an estimate of the information 
content of the satellite data at any given pressure 
altitude. The time series plots also reveal any sig- 
nificant biases between the data sets. Secondly, we 
present the monthly differences for each station in 
both tabular form and graphical form, but we do 
not address their statistical significance because of 
the small sample sizes. Near-seasonal statistics are 


generated at each station and compared with esti- 
mates of accuracy for each of the data sets. We also 
compare the 7-month statistics with the results in 
figure 11 of Gille et al. (1984a). 

The findings arc' discussed according to groups of 
pressure levels: (a) 10 hPa; (b) 5, 2. and 1 hPa; and 
(c) 0.4 and 0.1 hPa. We determine the statistical dif- 
ferences for the US versus USSR ROCOB’s and com- 
pare our findings with those of Gille et a.l. (1984b), 
who proposed using LIMS as a. transfer standard be- 
tween the US and USSR ROCOB's for 1978 79. To 
do this, plots of the monthly differences as a func- 
tion of station latitude are evaluated for the above 
pressure level groupings. Separate 7- mont h average 
statistics are presented for just the IFS comparisons. 

3. Data Sets 

3.1. LIMS Data 

The Nimbus 7 LIMS instrument was operational 
from October 25. 1978. to May 28, 1979. LIMS T(p) 
results on the LATPAT product were retrieved from 
C(> 2 radiances measured between 64° S and 84° N 
at approximately 4° latitude intervals, at a vertical 
resolution of about 2.5 km. but with a vertical point 
spacing of 1.5 km. The T(p) results were then inter- 
polated to 18 pressure levels from 100 to 0.05 hPaand 
synopticallv mapped to 1200 UTC using a Kalman 
filter technique. This Fourier coefficient product 
(termed LAM AT) was created at each standard lati- 
tude and pressure level (Remsberg et al. 1990). The 
standard NMC levels used here are a subset, of the 
LAMAT levels. 

The LAMAT data also contain geopotential 
height information Z(p) at. each pressure level, from 
which T(Z) or T(z) can be generated at a. station, 
where z is geometric altitude. The LIMS distribu- 
tions of Z(p) were obtained by use of the 50-hPa 
geopotential field as a reference and then integrat ing 
upward using the LIMS T(p) information. Because 
T(z) is a fundamental product of the Datasonde, one 
can also make LIMS/ROCOB comparisons of either 
T(z) or T(Z). 

Temperatures from the ascending and descending 
orbital segments (obtained at ~1 PM and ^11 PM 
at most latitudes) are different by less than 1 K in 
the lower st ratosphere. In order to have the full six 
zonal wave number or 30° longitudinal resolution in 
our LIMS analysis, we rely on the results obtained 
by combining data from all the orbital segments. 
Therefore, no provision has been made for any diur- 
nal temperature change. We note that diurnal varia- 
tions are significant in the upper stratosphere/lower 
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mesosphere at low latitudes, as estimated from sepa- 
rate zonal mean LIMS coefficients obtained at those 
two local times (Hitchman and Leovy 1985). Resolu- 
tion in the tangent-layer, limb-view direction is 200 
to 300 km or somewhat better than the 4° sam- 
pling resolution of the LAMAT data in the merid- 
ional direction. Remsberg et al. (1990) report that 
the LAIPAT temperatures have been mapped to an 
accuracy of about ±1 K. Thus, if there is no diur- 
nal variation or other bias in the original retrieved 
LIMS profiles, that value represents the average un- 
certainty of the mapped LIMS temperatures at a sta- 
tion location. 

3.2. In Situ Data 

In situ meteorological rocketsonde observations 
(ROCOB’s) provide high-resolution profiles of den- 
sity, temperature, and winds with altitude. Be- 
cause of the relatively high cost involved, however, 
ROCOB’s have been obtained routinely at only a few 
sites (e.g., about a maximum of 30 sites in 1965). 
Further, the frequency of observations at each sta- 
tion was about one sounding per week. We consider 
14 stations that were making soundings with Data- 
sondes during the LIMS period. The M-100 instru- 
ment was used at four USSR land-based sites. The 
set of 18 rocket stations used for the present study is 
listed in table 1. 

The Datasonde and M-100 temperatures are sub- 
ject to large errors at upper levels (Krumins and 
Lyons 1972; Nestler 1983). For example, major cor- 
rections were operationally applied to account for 
aerodynamic heating due to the rapid fall of each in- 
strument. Corrections for heat lag, radiation, and 
sensor emissivity were also included. For the US 
Datasonde these combined corrections are about 2 K 
at 40 km and 8 K at 60 km, while for the USSR 
M-100 they are much larger, particularly in the meso- 
sphere, as reported from a 1973 rocketsonde inter- 
comparison campaign (Finger et al. 1975). Because a 
large correction is less accurate, improvements were 
made to the M-100 payload (the M-100B system) 
for which the corrections are smaller, and a second 
intercomparison was conducted with the Datasonde 
in 1977 (Schmidlin et al. 1980). An IFS sensor system 
was also part of that 1977 intercomparison. Accord- 
ing to Schmidlin et al. (1980) and Koshelkov (1983), 
there was a gradual changeover by the USSR to the 
modified payload after 1978. But the printed copies 
of the profiles for our study still carry the M-100 sen- 
sor designation in the heading with no further com- 
ment on the nature of the corrections. Because we 
are not certain whether the 1978-79 USSR ROCOB 
temperatures were obtained with the original M-100 


or the improved M-100B payload design and their 
associated corrections, we consider comparisons with 
the USSR ROCOB’s to be qualitative at best. 

Magnetic tape versions of the ROCOB’s used in 
this study are available at the Wallops Flight Facility. 
Additional quality control criteria were not applied to 
each ROCOB, but, in general, the time series plots to 
be presented herein indicate that the use of strict ac- 
ceptance criteria would have been counterproductive. 
In fact, it is an inspection of the plots themselves 
that defines the quality of each of the data sets. A 
brief description of the T(p) profiles for each ROCOB 
type is given below. We have interpolated these T(p) 
profiles using cubic §pline techniques in logarithrq of 
pressure to give temperatures at the standard pres- 
sure levels for the present satellite validation study. 
The ROCOB profiles were interpolated rather than 
the satellite data because the vertical point spacing is 
no better than about 3.5 km for the mapped satellite 
data. In general, interpolation affects the random er- 
ror for a set of profile differences, but not its mean 
difference. 

3.2.1. US Datasonde. The US Datasonde in- 
strument, technique, and error sources are given in 
Schmidlin et al. (1980), Nestler (1983), and refer- 
ences therein. The precision or repeatability of the 
Datasonde T(z) is 1 K up to 53 km (Schmidlin 1981). 
However, above this level the Datasonde repeata- 
bility deteriorates exponentially to about 3.8 K at 
65 km and 7.5 K at 70 km. The pressure profile 
in the ROCOB was calculated by a tie-on of the 
rocket temperature-altitude T(z) profile to the geo- 
metric height derived from the 50-hPa level (nor- 
mally) of a colocated RAOB sounding, and then in- 
tegrated upward hydrostatically using the T(z) from 
the ROCOB. Occasionally, the rocket and RAOB 
profiles did not overlap in altitude, and in those cases 
there was an extrapolation of the RAOB data upward 
to achieve a tie-on point for pressure. In general, 
ROCOB’s from the US sites were made around local 
noon. 

Occasionally, a ROCOB T(p) is misregistered. 
One particular example is given in figure 1, which 
shows a LIMS/ROCOB comparison for May 7, 1979, 
at Thule. Note that the Datasonde T(p) in the 
left panel has its stratopause near 0.3 hPa. The 
right panel shows the Datasonde measurement in its 
more fundamental T(z) form. The colocated LIMS 
T(Z) profile is derived from the geopotential height 
field Z on the LAMAT product; thus there is a 
slight LIMS/ROCOB discrepancy because we have 
not distinguished between geopotential versus geo- 
metric altitude. Nevertheless, the temperature ver- 
sus “altitude” comparison is also poor. A colocated 


4 


RAOB sounding extends to 10 hPa and is nearly iso- 
thermal from there to 100 hPa, as is the LIMS 
profile. Deep isothermal layers are typical of the 
high-latitude lower stratosphere in spring. If the 
ROCOB/RAOB tie-on criterion is based on good 
temperature agreement at the lowest altitude of this 
ROCOB sounding (about 30 km), then that require- 
ment was met. In this case, because the RAOB 
sounding does not extend well above the 10-hPa level 
where the temperature is increasing, the observer 
cannot know that the ROCOB is incorrect, and not 
just anomalous. In fact, as one can see in appendix A, 
it is really the temperature time series comparisons in 
the upper stratosphere from Heiss Island and Thule 
that strongly suggest an altitude registration prob- 
lem for the ROCOB in figure 1. The problem is most 
likely due to uncertainties in angular pointing for the 
GMD-4 tracking system used at Thule at that time. 
(The more precise FPS-16 system was used at most 
US stations.) Misregistration becomes much less of 
a problem at low and mid latitudes, where the tem- 
perature profile varies with height even in the low 
to mid stratosphere. In those cases, misregistration 
would give a clear mismatch with a colocated RAOB, 
causing the ROCOB to be unusable or rejected. 

3.2.2. USSR M-100 . A brief description of the 
M-100 system and instrument is given by Schmidlin 
et al. (1980) and Finger et al. (1975). The M-100 
often had a measurement of static pressure from a 
Pirani heat manometer. The procedure for obtain- 
ing the final T(p) in those soundings involved iter- 
ating between both this measured pressure and a 
hydrostatic calculation of the pressure profile until 
a match was achieved. Schmidlin et al. (1980) dis- 
cusses two data- processing methods considered by 
the USSR. The “standard' 7 processing method was 
used for the M-100 ROCOB’s and is based on samples 
obtained every 30 sec during descent, which means 
that the corrections applied to them operationally are 
not very accurate in the mesosphere. The “prospec- 
tive” processing method was used with data obtained 
from the M-100B sensor design of the 1980 ’s, and it 
is based on samples taken every 5 sec. 

We decided to include the temperatures from the 
four USSR stations in this report because the indi- 
vidual time series plots are informative and because 
we can use them to assess any statistical differences 
between the overall set of US and USSR ROCOB’s. 
Generally, because of the high northern latitude loca- 
tions of two of the USSR sites, observation times fall 
during local nighttime there from November to the 
first week of February but change to local daytime 
from March to May. 


3.2.3. US Super Loki Sphere or IFS. As 
part of the ongoing intercomparison studies during 
1978-79, soundings with the IFS sensor were ob- 
tained at Ascension, Barking Sands, Wallops Island, 
Cape Canaveral, W r hite Sands, and Point Mugu. 
These occasions are identified by an “x” on the 
points in the time series plots. Those soundings can 
be used to test the adequacy of the corrections being 
applied to the Datasonde temperatures in the meso- 
sphere. The IFS p(z ) profile is obtained from the den- 
sity profile by integrating the hydrostatic equation 
downward from an assumed state at the top (near 
90 km), and then temperature is obtained using the 
gas law. Therefore, IFS T(p) is not dependent on a 
colocated RAOB sounding (Schmidlin 1984). 

After deployment for a given launch, the sphere 
is inflated to a superpressure of nominally 10 hPa. 
If complete inflation is not achieved, then the sphere 
becomes compressed prematurely at its lowest alti- 
tudes, causing the density profile (and inferred tem- 
perature) to be less accurate. 

3.3. National Meteorological Center 

(NMC) Data 

The NMC temperature data used here refer to 
analyses at 1200 UTC at stratospheric levels (Finger 
et al. 1965; NMC Office Note 84, “Packing and Iden- 
tification of NMC Grid Point Data,” June 1989). At 
10 hPa the 1978-79 NMC Northern Hemisphere anal- 
yses used RAOB data as input, but only after the 
data had been corrected for solar heating and radia- 
tive cooling effects (see below). Both the original 
and the corrected RAOB data are stored at the Na- 
tional Climatic Data Center (NCDC). Many of the 
original RAOB 7 s contain corrections applied at the 
stations or “pretransmission” corrections (Mclnturff 
et al. 1979; Mclnturff and Finger 1968). In general, 
these corrections were applied to soundings from the 
1978-79 VAISALA, Kew, A-22, and RKZ sondes, but 
there is no information in the W 7 orld Meteorological 
Organization (WMO) transmission code to let NMC 
(or any other user) know for sure that the correc- 
tion was applied. Still, an examination by NMC of 
the 0000 UTC and 1200 UTC transmissions for even 
those corrected soundings reveals differences that at 
times are of the order of 1 to 2 K. Therefore, NMC ap- 
plied a solar heating correction to make the daytime 
data compatible with the nighttime data. A correc- 
tion for long- wave cooling was also made at 10 hPa, 
but the effect was assumed the same for all the sonde 
types. 

Because fewer RAOB’s ascend to 10 hPa, espe- 
cially in polar winter, the 10-hPa analyses also relied 
on 50-percent persistence plus a 50- percent upward 
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regression derived using climatological temperature 
data (Finger et al. 1965). The temperature per- 
sistence was based on the previous and following 
0000 UTC RAOB reports. Where RAOB data ex- 
isted, much weight was put on them for the analy- 
sis. The Southern Hemisphere analyses at 10 hPa 
were based on a combination of operational satellite 
and RAOB data. The information content for the 
1978 -79 NMC analyses equatorward of 20° N or S is 
based on an extrapolation from about 20° latitude 
(Randel 1987). 

For both hemispheres the NMC analyses from 5 
to 0.4 hPa are based on satellite soundings. From 
September 24, 1978, to February 23, 1979, the NMC 
temperature analyses for 5 to 0.4 hPa were derived 
using regression equations (Gelman and Nagatani 
1977) based on radiances from the Vertical Tem- 
perature Profile Radiometer (VTPR) flown on the 
NOAA 5 satellite, and from February 25, 1979, to 
January 20, 1980, on radiances from the Strato- 
spheric Sounding Unit (SSU) on NOAA 6 (chan- 
nels 25 and 26). Both the VTPR and the SSU are 
nadir-viewing instruments and have a vertical res- 
olution of the order of 10 to 17 km. The weight- 
ing functions for SSU channels 25 and 26 have a 
vertical width of about 15 km and arc centered at 
about 15 and 6 hPa, respectively. Thus, the anal- 
yses are based on radiances from deep atmospheric 
layers. The regression relationships also depend on a 
climatological set of colocated rocketsonde/satellite 
soundings (Gelman et al. 1982; Anon. 1978; Gelman 
and Nagatani 1977). Because the USSR rocketsonde 
data in that climatology were made warmer in the 
mesosphere to make them compatible with the Data- 
sonde measurements, the 1978-79 NMC temperature 
analyses at 0.4 hPa are also dependent on Datasonde 
accuracy (Anon. 1978). Finally, as at 10 hPa, the 
NMC temperatures were extrapolated equatorward 
from 20° N or S at these higher analysis levels. 

In our study the NMC fields have been modified 
for easier use in the analysis of stratospheric tem- 
perature fields. The original gridded NMC data are 
represented by the coefficients of a harmonic series 
at 45 latitudes with a separation of 4° in latitude 
(88° S to 88° N). The stratospheric data sets were 
fit with 25 zonal coefficients (a zonal mean term plus 
12 sine and 12 cosine terms or 12 wave numbers), 
which gives a longitudinal resolution of 15°. Tem- 
peratures were then obtained at the exact longitude 
of each station for our comparisons, which, in effect, 
is an interpolation from the original 25 coefficients. 
Thus, the zonal resolution for the NMC temperatures 
is potentially better than that from the six zonal wave 
number LAMAT data, particularly at low and mid 


latitudes. Of course, our 1978-79 10 hPa compar- 
isons still depend on operational NMC analyses that 
are based solely on the number of RAOB reports 
available when the analysis was made. 

Gelman et al. (1982 and 1986) and Finger et al. 
(1993) compared NMC temperatures with ROCOB’s 
for periods when there was an operational changeover 
between two satellites with similar instruments. 
Based on those comparisons, they report tempera- 
ture adjustments for each of the pressure levels as 
a function of latitude. Thus, after adjustment it is 
easier to evaluate small trends in temperature over 
a period when a succession of two or more satellites 
was operating. For the present study we made the 
NMC temperature fields compatible by applying the 
recommended adjustments from 5 hPa to 0.4 hPa. 

4. Satellite/ROCOB Station Results 

The figures in appendix A are plots of LIMS, 
ROCOB, and NMC time series temperature compar- 
isons at the 10-, 5-, 2-, 1-, 0.4-, and 0.1-hPa pressure 
levels for all stations listed in table 1. Summary fig- 
ures in appendix A are also shown of the monthly 
mean T(p ) differences for LIMS minus ROCOB at 
each station. The monthly and 7-month summaries 
at each station are useful in assessing whether there 
might be a LIMS temperature bias that varies with 
pressure-altitude or station location. 

Statistics of temperature differences for LIMS mi- 
nus ROCOB, LIMS minus NMC, and NMC minus 
ROCOB were calculated for each station for the 
VTPR/LIMS (Nov.-Feb.) period, the SSU/LIMS 
(Feb. May) period, and the entire 7-month period. 
The monthly mean differences are also calculated 
for individual stations. November -February and 
February-May represent approximate winter and 
spring seasons, respectively, for the Northern Hemi- 
sphere. The winter period is characterized by strato- 
spheric warming activity, while the springtime at- 
mosphere is less perturbed and relaxing toward a 
radiative equilibrium state. Those “seasonal’' and 
7-month average results are tabulated in appendix B 
for each of the six (five for NMC) pressure levels. 

Previous satellite validation studies found some 
rather large differences with ROCOB measurements, 
especially at upper levels (for example, Gille et al. 
1984a and 1984b; Petzoldt 1979). We wanted to in- 
vestigate those differences in more detail. Results in 
appendix A also indicate some persistent, large bi- 
ases even at 10 hPa, particularly for the high-latitude 
USSR stations. Therefore, we applied a rejection cri- 
terion of temperature difference greater than ±20 K 
in order to edit out profile pairs that contain a poten- 
tially spurious result. At 0.4 and 0.1 hPa it could be 
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argued that even greater differences are likely, given 
the uncertainties for the several rocket measurement 
systems. 

In the “samples” column for the summary fig- 
ures in appendix A and the tables in appendix B, the 
quantity on the left side of the slash (/) is the num- 
ber of observations that entered into the calculations, 
while the quantity on the right is the number of ob- 
servations that exceeded our 20 K criterion. For the 
NMC/ROCOB comparison, the right-hand number 
also reflects any missing days in the NMG analysis. 
The 20 K criterion was not applied in LIMS/NMC 
comparisons. 

The mean differences for the paired observations 
and their standard deviat ions have been compiled in 
appendix B at each station and for each of the three 
time periods (7 months, Nov. Feb., and Feb. May). 
It is assumed that both observations (satellite and 
RQCOB) are representative of the same volume of 
the atmosphere, and it is our expectation that both 
observation techniques will register the true atmo- 
spheric temperature. Therefore, we have calculated 
quantities that can be used to test the hypothesis 
that the sample mean difference is zero (Guenther 
1973). The sample mean difference (or d) at. a sta- 
tion is given by 

Dj d i 


d = 


i) 


where 4 dj — Xj\ — xyi\ Xj\ and x yi are the satellite 
and ROCOB values, respectively, for pair j\ and rn is 
the number of pairs for that time period and pressure 
level. The corresponding standard deviations of the 
differences is defined from 


E Mi - df 

( rn - 1) 


( 2 ) 


The Student's t-test. statistic is defined as 

d 

t in - 1 - 


( 3 ) 


and can be determined from the quantities in 
appendix B. The quantity in the denominator of 
equation (3) is the standard deviation of the mean 
difference. For a given confidence interval, one can 
estimate whether the calculated d is significantly dif- 
ferent than zero, and thus whether there is a signif- 
icant seasonal bias in T(p) for at least one of the 
measurement techniques. 

The LIMS and NMC temperatures are for 
1200 UTC. For the US sonde stations, the obser- 
vations are taken near local noon. However, the 


time difference with ROCOB’s can be as much as 
12 hr, depending on the longitude of a station (e.g., 
Kwajalein -8.7° N, 168° E). The observations for 
USSR sondes at Volgograd and Heiss Island are 
in darkness from November to the first, week of 
March. The Thumba observations are within 4 hr 
of 1200 UTC. 

4.1. 10 hPa 

In general, the roeketsonde measurement errors at 
10 hPa are small and the differences from the LIMS 
and NMC analyses are small (see tables in appen- 
dix B). In the tropics the temporal small-scale vari- 
ability is more pronounced in the ROCOB’s than in 
the LIMS or NMC results. Some of the ROCOB 
variability may be due to the tie-on uncertainties for 
Datasonde T(p) profiles. Also the ROCOB s con- 
tain effects of small-scale oscillations due to tides and 
gravity waves. The damped amplitudes in the LIMS 
and NMC temperature time series in the figures for 
10 hPa in appendix A are attributed primarily to 
their lower zonal resolutions and to constraints in 
the map analysis products. At Molodezhnaya and 
Heiss Island and to a smaller extent at Fort Sher- 
man, Kwajalein. Thule, Thumba, and Volgograd, 
LIMS is almost always warmer than ROCOB’s. At 
Wallops Island LIMS is colder than ROCOB’s for 
most months. Four of those eight stations obtained 
ROCOB’s with the M-100 system for which signifi- 
cant corrections were applied routinely at the sites, 
particularly at the upper levels. Part of the compari- 
son biases may also be due to the finite (2.5 km) verti- 
cal resolution for LIMS, but. if so, it should be notice- 
able for all low- and mid-latitude stations, regardless 
of season. We have obtained LIMS retrievals with a 
higher resolution algorithm, and they do give tem- 
peratures that an 4 colder by about 0.5 K to 1.0 K at 
10 hPa for a sample day, January 13, 1979 (Solomon 
ct al. 1986). 

A slight bias in ROCOB temperatures (or, more 
likely, its pressure registration) would show up most 
clearly in tropical T(p) data because of their strong 
vertical gradients in the stratosphere. Diurnal tem- 
perature variations can also be a factor at low lat- 
itudes. Hitchman and Leovy (1985) found day 
temperatures colder than night by up to 1.4 K in 
the zonally and 2 16- day- aver aged LIMS results near 
10° N and 10 hPa: this difference is related to the 
semiannual oscillation and is most pronounced in 
Northern Hemisphere spring. LIMS temperatures 
used in the present study are merely an average 4 of 
the local day and night values. Consider the appar- 
ent biases at Kwajalein (fig. A4) and Fort Sherman 
(fig. A5), where ROCOB s were taken near midday. 
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local time. The 1200 UTC LIMS results in those fig- 
ures ought to be too warm due to not accounting for 
this diurnal temperature tide. The comparisons at 
Thumba (fig. A3) may be affected also. 

The NMC time series comparisons with ROCOB's 
are similar to those for LIMS at the individual sta- 
tions in appendix A, except that the NMC tempera- 
tures show less seasonal variability equator ward of 
about 10° latitude. The NMC analyses for Fort 
Sherman in figure A5 are a bit colder than the 
ROCOB’s — in opposition to the corresponding LIMS 
result. Of course, we are mindful that the NMC 
T(p) analyses are based on RAOB data that include 
bias “corrections” for solar radiative heating effects 
at low sun angle (sunrise), leading to a 1200 UTC 
NMC temperature at Fort Sherman that may be too 
cold. The 10-hPa NMC results for Thumba contain 
almost no “short” period variations, most likely be- 
cause very few RAOB reports from that region of the 
world were incorporated into the operational NMC 
analyses (Randel 1987). 

For mid- and high-latitude stations both LIMS 
and NMC contain the large-scale temperature vari- 
ations also seen in the ROCOBs. LIMS and NMC 
faithfully reproduce the warming events of December 
and January at Fort Churchill, Poker Flat, Primrose 
Lake, and Volgograd. These findings for LIMS agree 
with those from the LIMS/RAOB time series com- 
parisons for Berlin (52° N) at 10 hPa in Grose et al. 
(1988). Some biases remain for the polar stations. 
For example, Gille et al. (1984b) reported LIMS 
warmer than ROCOB’s at Heiss Island by about 6 K 
on average and warmer than ROCOB’s at Thule by 
about 2 K. The comparisons in appendices A and B 
are in accord with their findings. Even so, the qual- 
ity of the high-latitude LIMS results is judged bet- 
ter than for NMC, because the standard deviations 
for the LIMS/ROCOB differences are almost always 
smaller than those for NMC/ROCOB. 

The monthly mean differences for LIMS minus 
ROCOB’s are plotted as a function of station latitude 
in figure 2 for 10 hPa. The three Northern Hemi- 
sphere USSR stations are marked by open circles 
(“O” )• Sample size per month is small for all sta- 
tions and is given to the left of the “slash” at the right 
margins of each plot. Those samples were included in 
the final statistics; samples to the right of the “slash” 
were rejected. The last panel in the sequence in fig- 
ure 2 is the 7- month statistics, where the horizontal 
bars represent the standard deviations with respect 
to the MD for the 7-month period. There is no clear 
latitudinal trend in the LIMS/ROCOB comparisons. 


We focus on those instances where persistent 
monthly biases can be noted from the summary 
figures in appendix A; the 10-hPa LIMS/ROCOB 
statistics at individual stations are given in appen- 
dix B for two seasonal periods. For example, ta- 
ble Bll for Wallops Island has a value of d of —4.2 K 
for the November-February period with s = 3.1 K 
and m - 13 (see eqs. (1_) and (2)). Table B12 for 
Volgograd has a value of d of 2.5 K for the February- 
May period with s = 2.0K and m = 40. Table B5 for 
Fort Sherman has a value of d of 2.4 K with s = 2.5 K 
and for m = 32 for the 7-month period. Finally, ta- 
ble B18 for Heiss Island has a 7-month value of d of 
5.4 K with s — 3.4 K and m = 65. In each case, the 
differences, according to equation (3), are significant 
at the 99-percent confidence level. 

Figure 3 shows the NMC/ROCOB results, and 
the MD’s are similar in magnitude but opposite in 
sign at low latitudes to those for LIMS/ROCOB 
in figure 2. The 7-month average standard devia- 
tions are about equal. There are significant biases at 
Molodezhnaya (table Bl, Feb. -May), Thumba (ta- 
ble B3), Wallops Island (table Bll, Nov.- Feb.), and 
Heiss Island (table B18, Feb.-May). NMC/ROCOB 
standard deviations are larger than LIMS/ROCOB 
values at Fort Churchill, Thule, and Heiss Island. 

Figure 4 shows LIMS/NMC differences at 10 hPa. 
Sample size is much greater here, comprising essen- 
tially all days of each month. There are pronounced 
and persistent differences at low latitudes. Compar- 
isons with ROCOB’s are closer for LIMS at Ascen- 
sion, Kwajalein, and Thumba, but closer for NMC at 
Fort Sherman. Randel (1987) notes that the 1978-79 
NMC analyses equatorward of about 20° N or S are 
based on an extrapolation from 20° latitude, a pro- 
cess that may be less accurate over Asia where almost 
no RAOB reports were available for the operational 
analyses. 

At Ascension (8° S) the LIMS/ROCOB and the 
NMC/ROCOB differences are not significant. This 
finding is at odds with that for Ascension and for 
Natal, Brazil (6° S, 325° E) in Barnes et al. (1991), 
where they found NMC warmer than ROCOB by 
about 6 K and 7 K, respectively, at 10 hPa in spring 
1985. We note that NMC no longer made use of 
RAOB data in their 10 hPa Southern Hemisphere 
analyses after October 16, 1980, but relied solely 
on TOVS analyses (Gelman et al. 1986). More im- 
portantly, March 1985 was a transition period be- 
tween NOAA 7 and NOAA 9 for the NMC analyses. 
Although the so-called NMC adjustment factors at 
5 hPa were different by 3.8 K for those two satellites, 
no factors were developed for 10 hPa even though 
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temperatures at that level are based on TOVS SSU 
data, too. 

The LIMS/NMC differences at Thule (77° N) and 
Heiss Island (81° N) in January and February (fig. 4) 
are relatively large, but there is good agreement 
between them in spring. Even so, both the LIMS 
and NMC comparisons with the M-100 at Heiss 
Island show differences that are clearly positive in 
spring (figs. 2 and 3). LIMS and NMC comparisons 
with Datasondes at Thule are not as consistent, 
particularly for November-February. When NMC is 
compared with Thule ROCOB’s (table B17), there 
is no significant seasonal bias, although the SD is 
quite large. As noted in the discussion of figure 1, 
the altitude registration is also not accurate for some 
springtime Datasonde soundings at Thule. 

4.2. 5, 2, and 1 hPa 

The LIMS minus Datasonde time series compar- 
isons in appendix A at 5, 2, and 1 hPa are similar 
in value and character to the comparisons at 10 hPa 
(see summaries in figs. 5, 6, and 7). There is no 
clear bias with latitude. Station standard deviations 
in appendix B are larger near the stratopause be- 
cause of increasing T(p) measurement and colocation 
uncertainties as well as effects from tides and grav- 
ity waves. There are significant biases in winter at 
5 hPa for Barking Sands (3.3 K) and Wallops Island 
(—4.4 K) and at 2 hPa for Primrose Lake (4.6 K). 
In spring there is a bias at 5 hPa for Primrose Lake 
(2.9 K) and at 2 hPa for Fort Churchill (5.5 K). But 
there is no springtime bias at Shemya or Poker Flat 
at 2 hPa. An inspection of individual profile com- 
parisons at Fort Churchill reveals a sharp decrease in 
the Datasonde T(p) values from 1 hPa to 2 hPa that 
is not followed so well by LIMS because of its finite 
vertical resolution. 

The LIMS minus M-100 results have larger biases, 
and they stand out in the 7-month summary plots in 
figures 5, 6, and 7. Measurements at Thumba (ta- 
ble B3) show a 7-month difference that increases from 
0.3 K at 5 hPa to 4.4 K at 1 hPa. Part of this dif- 
ference profile can be explained by diurnal tempera- 
ture variations (Gille et al. 1984b). Measurements at 
Thumba station were taken near twilight or at night 
(local time of 7 PM to 12 PM). Differences for the two 
high-latitude stations are significant and consistently 
positive for both seasons. This is a clear indication 
of bias, most likely due to an overcorrection for the 
large aerodynamic heating term in the reduction of 
the M-100 sensor data. The November-February bias 
of about 10 K at 1 hPa for Volgograd (table B12) 
and Heiss Island (table B18) agrees closely with the 


recommended 8 K adjustment for the correspond- 
ing winter period (May- Aug.) for M-100 data at 
Molodezhnaya at 68° S (see table 3 in Koshelkov 
1983). 

NMC minus ROC OB comparisons are summa- 
rized in figures 8, 9, and 10 and in appendix B. 
NMC temperatures at these levels are based on satel- 
lite data (VTPR or SSU). In contrast to LIMS mi- 
nus ROCOB comparisons, monthly differences at sta- 
tions are variable and have a tendency to change sign 
with altitude (and perhaps atmospheric state). Sta- 
tion standard deviations are largest during the win- 
ter months. Nadir satellite temperature sounders 
have vertical resolutions of the order of 17 km in 
the upper stratosphere (Peckham 1974; Nash and 
Forrester 1986; Jackson et al. 1990). Under disturbed 
atmospheric conditions, NMC’s use of regression of 
the observed VTPR (or SSU) radiances against a 
climatology of rocket profiles may misrepresent at- 
mospheric temperature at a given pressure-altitude, 
even though the deep-layer averaged temperatures 
are accurate (Gelman and Nagatani 1977). 

An interesting example of this insensitivity to 
real atmospheric variations occurs at White Sands 
during the second half of December 1978. The 
5-hPa time series plot (fig. A9) shows a cooling 
trend from December 12 to 17, reaching a minimum 
before starting to warm up again, according to both 
LIMS and ROCOB’s. A nearly opposite trend is 
seen at 2 and 1 hPa in both LIMS and ROCOB’s 
during the same time period. However, the finer 
structure present in the higher resolution LIMS and 
rocket data is absent in the NMC analysis. NMC 
temperatures at all three levels are nearly constant 
during December at White Sands; NMC is about 
25 K colder than ROCOB J s in mid December at 
1 hPa! In fact, the NMC statistics that we report at 
1 hPa are actually better than they should be because 
3 of the 17 measurements in December exceeded 
our 20 K cutoff criterion and were not included 
in the seasonal difference. Similar problems for 
the December 1978 NMC analysis occur at four 
other mid-latitude stations (Barking Sands, Cape 
Canaveral, Point Mugu, and Wallops Island). 

NMC temperatures for this December period 
were derived from VTPR channels 1 and 2. Chan- 
nel 1 measurements are centered at 30 hPa, and 
80 percent of its energy conies from the 100- to 2-hPa 
region. Channel 2 peaks at 10 hPa, and 80 percent 
of its energy comes from the 100- to 5-hPa region 
(Gelman and Nagatani 1977). This smearing of en- 
ergy over such a wide altitude range can lead to in- 
accurate analyses at 1 hPa. In another example, 
there is also a substantial NMC/ROCOB bias at 
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Poker Flat at 1 hPa for late January. However, the 
NMC VTPR analyses tend to follow the ROCOB 
temperature time series at high latitudes better, most 
likely because variations in satellite nadir radiance 
measurements are of larger amplitude and occur over 
deep layers for winter at high latitudes. 

Over the SSU period (or spring 1979 for these 
comparisons), NMC temperatures were obtained 
from channels 25 and 26. It must be stressed that 
channel 27 centered near 1.5 hPa— -was not oper- 
ational for this particular SSU instrument, caus- 
ing some degradation in NMC temperature accu- 
racy at 2 and 1 hPa. For example, this may be 
the reason for the large and statistically significant 
NMC/ROCOB bias for Poker Flat at 2 and 1 hPa 
during March, April, and May, when wave activ- 
ity was weak. NMC displays excellent agreement 
at 2 hPa with the springtime ROCOB’s at Fort 
Churchill, but not at 1 hPa. Most likely, the clima- 
tological profiles, used for regression by NMC, have 
a shape that is also different from the real atmo- 
sphere at that station (as defined by the ROCOB's 
for 1979). The transition from VTPR to SSU in the 
NMC analyses occurs on February 23, but any re- 
maining uncorrected discontinuities in the NMC time 
series are hard to distinguish from the temperature 
fluctuations that also occurred then. 

Another interesting result can be seen at 1 hPa 
for the higher latitudes in the panels for March and 
April. Figure 7 shows that L1MS matches the Data- 
sondc very well near 50° N and 80° N, but not 
M-100. Conversely, figure 30 shows that NMC tends 
to be colder than Datasonde but warmer than M-100. 
Because the NMC results are constrained more by 
a ROCOB climatology at this time, it is reason- 
able that the NMC analyses reflect that climatology. 
But since the high-latitude M-100 results at 1 hPa 
are significantly colder than LIMS (by about 10 K) 
and since NMC applied an adjustment of only 2 K 
to the M-100 data at 50 km when they compiled 
their rocket climatology (Anon. 1978), it is likely 
that their high-latitude climatology defines an atmo- 
spheric state based on both the Datasonde soundings 
and the undercorrected M-100 soundings. Hence, the 
retrieved NMC T(p) values at 1 hPa ought to split 
the differences between the two rocket sensors, as it 
seems to do in figure 10. 

The results of NMC minus Datasonde compar- 
isons corresponding to the two different satellite 
periods are reported in Gelrnan et al. (1982). Gen- 
erally, the station MD’s in appendix B are smaller 
at 5 and 2 hPa for both periods than in that refer- 
ence. This is to be expected because, as stated earlier 
in section 3, our NMC results do incorporate their 


recommended adjustment factors to make the differ- 
ent NO A A satellite measurements compatible. How- 
ever, our SD values are still comparable with those 
from Gelrnan et al. (1982) for both the corresponding 
VTPR period and the SSU period. 

Finally, a summary of the LIMS/NMC compar- 
isons is provided in figures 11, 12, and >3. In gen- 
eral, there is better agreement at low latitudes at 5 
and 2 hPa than at 10 hPa (compare fig. 4). However, 
a bias appears for the mid-latitude stations, leading 
to a distinct latitudinal dependence in the monthly 
plots. The character of that bias at 2 hPa is very 
similar for each of the SSU months of March to May. 
The station SD values in appendix B are small, in 
spring compared with winter, when the atmosphere 
is more variable. Several locations have statistically 
significant biases. In particular, there is a bias for 
the four stations from 29° N to 38° N for March at 
5 hPa (fig. 11), which should not be related to the 
loss of the top SSU channel. For comparison, figure 5 
shows only weak LIMS/ROCOB biases at Wallops 
Island and Point Mugu and none at White Sands. 
Figure 8 shows larger negative NMC/ROCOB bi- 
ases at all three stations. It is likely that this ap- 
parent NMC discrepancy in figure 11 is related to 
the coarse vertical resolution of the SSU and the re- 
gression procedure used to derive T(p) values from 
its radiances. Those constraints could also account 
for a high-latitude LIMS/NMC springtime bias at 2 
and 1 hPa (figs. 12 and 13). 

4.3. 0.4 and 0.1 hPa 

LIMS/ROCOB comparisons were made at 
both 0.4 and 0.1 hPa (appendices A and B), and the 
summary results are provided in figures 14 and 15. 
The 7-month LIMS minus Datasonde comparisons 
show that LIMS is a bit cold at 0.4 hPa, but clearly 
so at 0.1 hPa. Figures 14 and 15 indicate no clear 
seasonal or latitudinal dependence in those differ- 
ences. A top-of-profile effect may cause the LIMS 
temperatures at 0.1 hPa to be too cold by up to 2 K, 
but only when mid- mesosphere temperatures are cold 
enough to cause the radiance signal to approach the 
LIMS noise level. This retrieval bias occurs because 
an isothermal guess temperature is used initially at 
profile top. The effect of this guess disappears af- 
ter several iterations, except within about 3 km from 
the profile top. Retrievals generally begin at 0.05 
to 0.08 hPa. At any rate, this effect does not explain 
the large negative bias (—9 K) at 0.1 hPa in figure 15. 

LIMS temperatures are warmer than those from 
the M-100, an effect that is the opposite of that with 
the Datasonde. Mean differences with the M-100 
decrease progressively from about 8 K at 1 hPa to 
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about 4 K at 0.1 hPa. The SD increases, however, 
from 1 hPa to 0.1 hPa. The 7-inonth summaries 
in figures 14 and 15 show that both the sample 
differences and the standard deviations increase with 
latitude at 0.4 and 0.1 hPa. 

The highest analysis level for NMC is 0.4 hPa. 
Time series plots for each station (in appendix A) 
show a nearly constant NMC result at that level. 
There is also no SAO signature at low latitudes. 
On the other hand, there is no bias in the 7-month 
results for NMC minus Datasonde (fig. 16). This 
may be because the VTPR and SSU channels are in- 
sensitive to atmospheric temperatures from that 
level, such that the NMC regression procedure relies 
almost entirely on the long-term Datasonde climatol- 
ogy. There is a relatively low correlation coefficient 
(0.55) for the VTPR/ROCOB regression relation for 
temperature at 0.4 hPa (Anon. 1978). Presumably, 
that coefficient would be even smaller for the present 
SSU period with its top channel missing. 

The transition from VTPR to SSU in late Febru- 
ary in the NMC analyses is marked by a noticeable 
increase in temperature at Thule and Heiss Island 
(appendix A) at 1.0 and 0.4 hPa. This increase 
is opposite the direction of the temperature trend 
recorded by both ROCOB’s and LIMS. After the 
nearly 2-week transition period, NMC does follow the 
ROCOB’s more closely. 

Figure 17 indicates pronounced LIMS minus 
NMC differences at 0.4 hPa both by month and by 
latitude. Based on the rocket comparisons, one might 
conclude that both LIMS and NMC have significant 
errors, but of course, that conclusion also depends on 
the accuracy of the rocket T(p) in the mesosphere. 
Nevertheless, the individual station rocket /satellite 
time series analyses in appendix A do show that the 
rocket temperature trends were followed much better 
by LIMS than by NMC, indicating that the LIMS 
temperatures are more precise. 

4.4. Comparisons at 68° S 

Appendices A and B also include results for the 
USSR station Molodezhnaya (68° S, 46° E). Because 
the LIMS data do not extend beyond 64° S, the 
LIMS/ROCOB comparisons are not as useful for val- 
idation purposes, especially in autumn (April and 
May) when the Southern Hemisphere polar vortex is 
well formed and there is a larger meridional temper- 
ature gradient at high southern latitudes. To see this 
better, we have included two NMC curves one for 
64° S and another for 68° S. While the effect of the 
gradient is apparent, there is still a significant bias 
with the rocket data in autumn. The NMC/ROCOB 


statistics in appendix B were calculated using NMC 
data at 68° S, and they show a large bias, too. 

During summer, when the gradients in the mid 
stratosphere are weaker, there is very good agree- 
ment between LIMS, NMC, and the rocket data at 
5 and 10 hPa. Both LIMS and NMC are warmer 
than the M-100 at 2, 1. and 0.4 hPa, as was the case 
for the other three USSR rocket stations. The LIMS 
T(p) time series at 0.1 hPa is nearly constant from 
October to April with a gradual warm up in autumn. 

5. Comparisons by Latitude Zone 

5.1. LIMS/Datasonde 

This section focuses on the LIMS/Datasonde 
comparisons. We have grouped those paired results 
by latitude zone to search further for any mean differ- 
ences. There an' five stations at low latitudes (8° S 
to 22° N). four at mid latitudes (29° N to 38° N), 
and five at high latitudes (53° N to 77° N). Equa- 
tions (1) and (2) have been applied to these larger 
samples, and the results are given in table 2 and fig- 
ures 18, 19, and 20. Horizontal bars represent the 
sample standard deviations from table 2. Mean dif- 
ferences are judged significant if they are outside the 
95-percent confidence intervals for t hese larger sam- 
ples. The important issues are how do the mean dif- 
ferences compare with the estimates of systematic 
error in T(p) for LIMS and Datasondes and do they 
vary with latitude', season, or pressure. 

LIMS is significantly warmer at 10 hPa by about 
1 K for low and high latitude, according to equa- 
tion (3). The high-latitude bias is similar for winter 
and spring, indicating that the problem is not likely 
due to a misregistration of any Datasonde profiles. A 
bias of 1 K is of the order of the expected accuracy 
of T(p) at 10 hPa for both LIMS and Datasonde. A 
LIMS bias of this order could be attributed simply 
to uncertainties in the transmit t a, rices for C(>2 (Gille 
et al. 1984a), although such a bias is expected to be 
fairly uniform with latitude. 

For the pressure range, 5 hPa to 1 hPa, there is 
a significant LIMS/Datasonde bias at high latitudes, 
but not at low or mid latitudes. The bias is most 
pronounced at 2 hPa in winter, when LIMS is warmer 
by 3.6 K (table 2). The standard deviation of that 
mean, s/y/rri, is 0.5 K, so the 95-percent confidence 
interval is about twice that, only ±1.0 K. On average, 
LIMS is warmer at all three pressure levels at high 
latitudes in both winter and spring. 

At 0.4 and 0.1 hPa LIMS is cooler than Datasonde 
for all three latitude zones. The mean difference 
of about -9 K at 0.1 hPa is also much greater 
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than the theoretical LIMS, root-sum-square (rss), 
T(p) error estimate of ±4.6 K in Gille et al. (1984a, 
their table 2). It is also larger than the remaining 
uncorrected Datasonde biases reported by Nestler 
(1983). 

Gille et al. (1984a) reported mean and standard 
deviation differences for three stations — Ascension, 
White Sands, and Fort Churchill. Our station results 
from appendix B are similar, although we do show 
smaller mean differences at W T hite Sands and larger 
ones at Fort Churchill. From the theoretical LIMS 
rss errors in Gille et al. (1984a), the T(p) error at 
10 to 1 hPa for a single profile is somewhat smaller 
than the 7-month SD values in our table 2. But then 
our SD values also include any uncertainties in the 
Datasonde T(p). Our 7-month MD values at low 
and mid latitudes arc much smaller than the LIMS 
rss errors, indicating that many of the systematic 
LIMS T(p) errors are quasi-random when averaged 
over many profiles. 

5.2. LIMS/Sphere 

The falling sphere (IFS) technique has been used 
for T(p) validation, primarily in the mesosphere. 
Schmidlin et al. (1991) indicate average IFS minus 
Datasonde differences in T(z) of less than 3 K from 30 
to 60 km with a repeatability of order ±3 K. Their 
results are shown in figure 21. The LIMS/Datasonde 
comparisons in figures 18 to 20 have shapes that are 
very similar to figure 21 and are just as accurate, 
even taking into account the additional estimates of 
LIMS T(p) error due to its finite vertical resolution 
and the fit of its mapped coefficients to the original 
retrieved profiles. 

Schmidlin et al. (1991) and Quiroz and Gelman 
(1976) did find Datasonde warmer than the IFS T(p) 
values by about 5 K at 0.1 hPa (about 65 km in 
fig. 21). It is believed that the IFS is the more accu- 
rate in situ technique in the mesosphere. Therefore, 
we use the Datasonde as the common data set for the 
two LIMS/in situ comparisons and find that LIMS 
minus IFS should be only about —4 K at 0.1 hPa. 
That difference is of the order of the rss error for the 
LIMS T(p). We infer then that the low- and mid- 
latitude LIMS T(p) values are accurate from 10 hPa 
to 0.4 hPa, but too cold at 0.1 hPa. 

As a check on our inferred LIMS/IFS differences, 
we calculated LIMS differences with the approxi- 
mately 70 IFS soundings that appear at 6 mid- 
latitude stations in the time series plots in appen- 
dix A. These IFS comparisons are more meaningful 
for us because they are collocated in time (same day 
and year), as well as space. The individual station 


and 6-station average results are given in table 3 
and figure 22 along with the 7-month mid-latitude 
LIMS/Datasonde plot from figure 19. Standard de- 
viations for LIMS/IFS are rather large, presumably 
because of the small-scale effects of the vertical winds 
in an individual sphere determination of T(p). Be- 
cause the standard deviation about the m$an is larger 
for the IFS comparisons, it is concluded that a T(p) 
profile from a single Datasonde may be more rep- 
resentative of the atmosphere than the T(p) from 
a single sphere. IFS sample size is smaller at 5 
and 10 hPa, because the sphere often deflates some- 
what before descending to those levels. Still, the 
LIMS/IFS differences are within ±3 K over most of 
the pressure range.' At 0.1 hPa LIMS minus IFS is 
equal to —4.6 K ± 4.8 K, a result that is remarkably 
similar to our inferred LIMS/IFS bias based on the 
findings in Schmidlin et al. (1991) and based on the 
LIMS/lidar comparisons in Remsberg (1986). 

6. Discussion and Conclusions 

The present time series approach to satellite tem- 
perature validation has enabled us to make use of all 
the rocketsonde data, thus increasing both sample 
size (to 665) and statistical confidence. Furthermore, 
one can more easily judge the quality of the measured 
temperatures at a station by observing the general 
agreement in their variations for both quiet and dis- 
turbed atmospheric conditions. The high precision of 
Limb Infrared Monitor of the Stratosphere (LIMS) 
temperature- versus-pressure profile (T(p)) values is 
particularly evident in our station time series plots. 

There are no clear trends with season for LIMS 
minus Datasonde T(p) for any latitude zone, indi- 
cating that the LIMS T(p) retrieval accuracy is not 
a function of atmospheric state. This is an impor- 
tant point, because it confirms the robustness of the 
LIMS temperature retrieval technique. There is no 
stratospheric LIMS bias at low and mid latitudes. 
However, we do find a positive LIMS bias for both 
seasons at high latitudes of the upper stratosphere, 
and it is hard to imagine how the Datasonde could 
be accurate at low and mid latitudes, but not at high 
latitudes. It is also unlikely that a latitudinally vary- 
ing LIMS bias could be explained by transmittance 
errors or a smoothing of the true T(p) due to the 
vertical resolution of LIMS. We have found a minor 
error (about 0.25 percent) in our first guess for the 
altitude above the center of the Earth of the tan- 
gent layer for our reference pressure level. Because 
this value is used in the hydrostatic calculation of 
the LIMS T(p), it leads to a systematic T(p) error 
that is increasingly positive at high latitudes of the 
upper stratosphere. The effect of this error for a 
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mesospheric T(p) is less clear and more difficult to 
sort out, at least by comparison with the Datasondes 
at high latitudes. (We have no sphere profiles for 
high latitudes during the LIMS period.) More work 
is needed with an improved LIMS algorithm plus up- 
dated spectral databases for the LIMS CO 2 channels 
to evaluate these effects in detail. 

Gille et al. (1984b) proposed using LIMS T(p) 
data as a transfer standard between the Datasonde 
and the M-100. We found differences in this study 
that are very similar to those of Gille et al. (1984b, 
their figs. 2, 3, and 4). Those differences also 
agree qualitatively with the recommended adjust- 
ments to the M-100 temperature climatologies for 
high-latitude stations, at least at 50 km (Koshelkov 
1983). However, the positive LIMS bias, noted 
above, would affect its use as a transfer standard at 
high latitudes. 

The 7-month National Meteorological Center 
(NMC)/Datasonde comparisons at 10 hPa show 
good agreement, except at low latitudes, where 
the 1978-79 NMC analyses are based on extrapo- 
lations of RAOB data equatorward of about 20° N. 
NMC/Datasonde mean differences are very similar 
to those for LIMS/Datasonde at 5 hPa, a level where 
nadir-radiance data were available to NMC during 
1978-79. SD values at 5 hPa, however, are larger 
for NMC as compared with LIMS, indicating that 
the LIMS analyses follow the true temperature vari- 
ations better than the NMC analyses. At 2, 1, and 
0.4 hPa the 1978-79 NMC retrievals are weighted 
toward their historical ROCOB climatology, and, as 
a result, the NMC/Datasonde time series compar- 
isons show larger differences at those levels. The 
largest differences occur during winter and are at- 
tributed to the low vertical resolution of the nadir- 
viewing sounders and a nonrepresentative climatol- 
ogy, as shown in the White Sands time series. Several 
stations (e.g., Poker Flat) also have relatively large 
NMC/ROCOB mean differences during springtime 
at 2 and 1 hPa, perhaps because of the relative in- 
sensitivity of SSU channels 25 and 26 to stratopause 
temperatures. 

The findings herein suggest that there is no sta- 
tistically significant T(p) bias affecting the LIMS 
species in the upper stratosphere at low and mid 
latitudes. The T(p) bias at high latitudes affects 
LIMS species there in two ways. First, there is a 
bias in the registration of the measured species ra- 
diance profiles with pressure. Second, a T(p) bias 
affects the calculation of blackbody radiances, which 
must be accounted for in a limb emission retrieval. 
Both effects, while small, must be evaluated further. 
Finally, because LIMS temperatures seem to be too 


cold at pressure levels above 0.4 hPa at all latitudes, 
the LIMS ozone values may be too large at those 
levels. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
January 5, 1994 
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Table 1. Rocketsonde Stations 


Serial 

no. 

Station 

identification 

Latitude, 

deg 

Longitude 

(east), 

deg 

Station 

name 

Instrument 

type 

1 

89542 

-67.7 

46 

Molodezhnaya 

USSR M-100 

2 

61902 

-8.0 

346 

Ascension Island 

US Datasonde 

3 

43373 

8.5 

77 

Thumba 

USSR M-100 

4 

91366 

8.7 

168 

Kwajalein 

US Datasonde 

5 

78801 

9.3 

280 

Fort Sherman 

US Datasonde 

6 

78861 

17.1 

298 

Antigua 

US Datasonde 

7 

91162 

22.0 

200 

Barking Sands 

US Datasonde 

8 

74794 

28.5 

279 

Cape Canaveral 

US Datasonde 

9 

72269 

32.4 

254 

White Sands 

US Datasonde 

10 

72391 

34.1 

241 

Point Mugu 

US Datasonde 

11 

72402 

37.8 

285 

Wallops Island 

US Datasonde 

12 

34560 

48.7 

44 

Volgograd 

USSR M-100 

13 

70414 

52.7 

174 

Shemya 

US Datasonde 

14 

71124 

54.8 

250 

Primrose Lake 

US Datasonde 

15 

71913 

58.7 

266 

Fort Churchill 

US Datasonde 

16 

70192 

65.0 

213 

Poker Flat 

US Datasonde 

17 

04202 

76.6 

291 

Thule 

US Datasonde 

18 

20046 

80.6 

58 

Heiss Island 

USSR M-100 
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Table 2. LIMS Minus Datasonde Statistics by Latitude Zone 
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Figure 1. Temperature profile comparisons at Thule (May 7, 1979). 
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Figure 2. LIMS minus rocketsonde temperature differences by month and station latitude at 10 hPa. Horizontal bars 
in last panel represent standard deviation about 7-month average difference. 
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Figure 2. Concluded. 
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Figure 3. Concluded. 
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Figure 4. Concluded. 
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;urc 5. LIMS minus rockctsonde temperature differences by month and station latitude at 5 hPa. Horizontal bars 
in last panel represent standard deviation about 7-month average difference. 
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Figure 7. LIMS minus rocketsonde temperature differences by month and station latitude at 1 hPa. Horizontal bars 
in last panel represent standard deviation about 7-month average difference. 
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Figure 7. Concluded. 










US Datasonde 
USSR M-100 



6ep ‘epnme“| 


Figure 8. Concluded. 
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Figure 9. Concluded. 
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Figure 10. NMC minus rocketsonde temperature differences by month and station latitude at 1 hPa. Horizontal bars 
in last panel represent standard deviation about 7-month average difference. 














US Datasonde 
USSR M-100 
Mean 




39 


Figure 11. Concluded. 
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Figure 12. LIMS minus NMC satellite temperature differences by month and station latitude at 2 hPa. Horizontal 
bars in last panel represent standard deviation about 7-month average difference. 
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Figure 12. Concluded. 
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;urc 13. Concluded. 
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Figure 14. Concluded. 
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Figure? 15. Concluded. 
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Figure 17. LIMS minus NMC satellite temperature differences by month and station latitude at 0.4 hPa. Horizontal 
bars in last panel represent standard deviation about 7-month average difference. 
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Figure 17. Concluded. 
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Figure 18. Seasonally averaged LIMS minus Datasonde (US) temperature differences for low latitudes. 
Horizontal bars represent standard deviation about average difference. 
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Figure 19. Seasonally averaged LIMS minus Datasonde (US) temperature differences for mid latitudes. 
Horizontal bars represent standard deviation about average difference. 
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Figure 20. Seasonally averaged LIMS minus Datasonde (US) temperature differences for high latitudes. 
Horizontal bars represent standard deviation about average difference. 
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Figure 22. LIMS minus ROCOB temperature differences at mid latitudes. 


54 



Appendix A 

Temperature Time Series Plots and 
Statistics Plots at Each Station 

This appendix contains temperature (K) time 
series (day) plots as well as monthly mean differences 
at each of the rocketsonde stations listed in table 1 
of this report. The time series plots are for the 10-, 
5-, 2-, 1-, 0.4-, and 0.1-hPa levels. The solid lines are 
LIMS data, and the dotted lines are NMC data for 
the latitude and longitude of the rocketsonde station. 
The filled circles represent rocketsonde data; a filled 
circle with an “x v on it indicates a falling sphere 
data point. In figure Al for Molodezhnaya. the LIMS 
data are from 64° S (LIMS data do not exist south 
of 64° S), and the dashed line represents NMC data, 
at 64° S for comparison with LIMS. The dotted line 
is NMC data at the station latitude (68° S). 


The monthly average profile plots that follow the 
time series plots for each station represent the mean 
differences for LIMS versus Datasonde or M-100. 
Rocketsonde observations were not available for all 
months at some stations. The 7-month average 
difference is also given, where the horizontal bars 
are the standard deviations about those differences. 
The first number in the parentheses on the right- 
hand border of each figure indicates the total number 
of samples that went into the calculation of the 
statistics. The second number indicates the number 
of observations that were excluded on the basis of a 
rejection criterion of 20 K (see text). The numerical 
values for individual stations, used in preparing these 
plots, are given in appendix B. 

The NMC data do not exist above the 0.4-hPa 
level (i.e., at 0.1 hPa in this study). 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure A2. Ascension Island (8° S, 346° E). 
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Figure A2. Continued. 
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Figure A3. Continued 
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Figure A4. Continued 
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Figure A4. Continued. 
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Figure A5. Fort Sherman (9° N, 280° E) 
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Figure A5. Continued. 
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Figure A6. Continued. 
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Figure A7. Barking Sands (22° N, 200° E). 
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Figure A8. Cape Canaveral (29° N, 279° E). 











Figure A8. Continued. 
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Figure A9. Continued. 
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Figure A 9. Continued. 























>1 ‘0jnjBjedLU8i 


>i ‘0jrnBj0duj0i 










(LIMS minus rocketsonde temperature) 



BcjLj ‘ajnssajd 


Bdi| ‘ajnssajd 


Figure A 10. Concluded. 
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Figure A12. Volgograd, USSR (49° N, 44° E). 







310 



>j ‘0jniEJ0diu0i 


y{ ‘0jniBJ0diu0i 


101 


Figure A 12. Continued. 
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Figure A13. Shemya (53° N, 174° E). 
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Figure A 13. Continued. 
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Figure A14. Primrose Lake (55° N, 250° E). 
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Figure A 14. Continued. 
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Figure A 14. Continued. 
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Figure A16. Continued. 







Figure A 16. Continued. 
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Figure A17. Thule (77° N, 291° E). 
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Figure A 17. Continued. 







Figure A 17. Continued. 
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Figure A18. Heiss Island, USSR (81° N, 58° E). 
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Figure A 18. Continued. 























Appendix B 

Tables of Temperature Statistics for 
Each Station 

This appendix contains tables of temperature 
statistics (mean and standard deviation of temper- 
ature differences in K) at each of the 18 rocket- 
sonde stations listed in table 1. Each table contains 
statistical comparisons between (a) LIMS minus 
rocketsonde, (b) LIMS minus NMC, and (c) NMC 
minus rocketsonde temperatures. The NMC data 
used here are from a Fourier coefficient product de- 
rived from the archived NMC gridded analyses. 

The '‘Total of samples” is the number of samples 
used in computing the statistics. Along with the 
“Number rejected,” the two values constitute the 


total number of observations at a given level for the 
corresponding period. 

The row labeled “7 month” represents statistics 
for November 1, 1978, to May 28, 1979. The row 
“Nov. -Feb.” represents the VTPR period (Nov. 1, 

1978, to Feb. 22, 1979) for NMC data. The row 
“Feb. -May” represents the SSU period (Feb. 25, 

1979, to May 28, 1979) for NMC data in this study. 
During the transition from VTPR to SSU (Feb. 22 
25, 1979) for NMC, no data exists above 10 hPa for 
the days of February 22, 23, 24, and 25, 1979, even 
though the time series plots are continuous. 

In table B1 for Molodezhnaya, the LIMS data are 
at 64° S and NMC? is at station coordinates. The 
7- month average LIMS minus rocketsonde results are 
plotted in appendix A. 
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Table Bl. Statistics for Molodezhnaya (08 
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Table B3. Statistics for Thumba (9° N. 77° E) 
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Table B 7 . Statistics for Barking Sands ( 22 ° N, 200 ° E) 
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Table Bll. Statistics for Wallops Island (38° N, 285° E) 



03 















G 


— 

— 



— 

— 



— 




CU 

rG 


iO 

cq 

o 


o 

O 

o 


o 

o 

o 



Oh 

JG 


C0 

CM 

LO 


o 

o 

o 


o 

o 

CP 


rH 


d 

LO 

d 


d 

d 

d 


d 

d 

d 





id 

td 

td 


d 

d 

d 


o 

d 

o 


O 















d 














<e 















G 













1 

Ph 

^=: 


o 

00 

h- 


tO 

LO 

p 


p 


p 



Oh 


p 

O 

p 


r-H 

o 

p 


CO 

CO 

CO 

G 

TT 


LO 

CO 

id 


d 

id 

CM 



00 

lO 


G 

■^r 


LO 

td 

d 


d 

rH 

cd 

CM 


to 

oc- 

d 

X 

O 















d 








































- 

G 














.. 

G 













o 



00 

CO 

p 


LO 

lo 

p 


p 


p 


G 

0 

Oh 


o 

LO 

p 


LO 

05 

p 


p 

p 

o 


’- a 


CO 

d 

cd 


LO 

d 

CM 


LO 

d 

CM 



jG 


o 3 

TH 

cd 


LO 

to 

CM 


LO 

td 

d 

G 

1—1 














£ 

r—i 













’> 















> 














<v 

o3 














G 

G 













TJ 

0-. 


O 

p 

P 


p 

p 

p 


O 

p 

00 


"O 

Oh 


oc 

rH 

p 


p 

05 

TP 


X 

oq 

h- 

T? 

jG 


LO 

LO 

d 


cd 

d 

r-H 


id 

d 

cd 


T? 

-C 


cd 

LO 

CM 


cd 

d 

rH 


cd 

d 

CM 

g 

CM 














Lh 

G 

CM 











"O 















T3 














G 

G 














G 

G 













G 

Ph 


p 

p 

p 


p 

P 

p 


p 

rH 

p 


G 

Oh 


p 

p 

t-> 


p 

p 

to 


o 

CO 

O 

cn 

Jh 


CO 

CM 

cd 


d 

d 

CM 


td 

d 

LO 


C/5 

JG 


CO 

Tji 

CM 


cd 

d 

H 


id 

h- 

cd 


lO 















lO 














G 















G 


V- 












Dh 


p 

p 

p 


p 

p 

p 


p 

p 

p 



H 

jG 


p 

O 


p 

TT 

p 


p 

CM 

'oc 


o 


CO 

CO 

cd 


d 

d 

d 


o 3 

cd 

^f 3 



o 


CM 

cd 

CM 


cd 

cd 

CM 


cd 

d 

CM 































o5 















G 
















CM 

T 1 



O 

o 

O 


o 

o 

o 



Ch 


o 

a 

0 


53 

r~~ 









LO 

h- 

00 


o 

o 

o' 


o 

o 

o 



gg 


d" 

CM 

05 


o' 




> 

o' 

S- 


J 














. 

1—1 


to 

CM 

CO 










o 













^ s 


d 



























W 















rejected a' 





























G 

Oh 


o 

LO 

o 

CO 

o 

CM 

-M 

U 

G 

T3 

o 

©^ 

o 

CO 

G 

"C 

G 

g 


CM^ 

CM 

o 

o 

G 

G 

G 

Oh 

JG 


r-H 

o 

c. 

o 

05 

CJ 

'O 

o 

TP 


o 

CO 

o 

0 

05 

CM 

LO 

r- 

LO 

tp 

d 

a 

1-0 

CM 



O 

M 

a 

OO 

2 

00 

u 

CM 

1—1 


£ 

o 

G 

G 

"o7 

p 

d 

0) 

to 

CM 

CO 

o 

a 

00 

o 

QO 

X 

o 

tO 


CO 














05 



T3 








CJ 




u 

a> 

n 

o3 

Ph 

G 

O 

w 

o 

o 

o 

G 

.2h 

O 

o 

o 

o 

(— 

to 

CO 

CO 

Tf 

t-i 

Q> 

G 

0H 

G 

p 

Tf 

CO 


H— > 

.S 

o 

o 

0 

o 

u 

o 

lO 

LO 


gG 


to 

CO 

CM 


T* 


CO 


05 

o 

C5 


G 




C5 






■X 




fi 

c 

CO 


CM 

ffl 

o 

00 

o 

00 


CM 



G 

g 

GG 

a; 

r- 

OC 

a> 



co 



h- 

p 


H 

^H 

o 






G 



Jh 

bC 

O 

G 

rH 

o 

LO 


CO 

■oo 

o 

■00 

.s 

LO 

rH 




o 




0 




g 





































o 








X 

o3 


o 

o 

O 

CJ 

o 


o 


LO^ 

CM 

CO 

bO 


G 

l-H 

CM 



u 








Ph 

X 






©^ 




X 





o 

o 

o 

y s 


TP 

lO 

V 

-G 

CM 

g 

G 

CO 

CO 

CO 

CO 

CM 

o 

S 

00 

© 

CO 

00 

O 

G 

T3 

CO 

- 

o 

CM 

£ 

(— 

"a 

CU 

JG 

CM 

d 

.£ 

H 

o 

to 

CM 

CM 

00 

CO 

a 

§ 

Tp 

00 

s 

CO 

oo 

4H 

CJ 

p 

T3 

\ 

LO 

oT 

Tp 

CO 

S3 

X 

t+H 






z 

X 




o 




o 

G 

X 

R 








0 




G 

Ph 

H 

o 

o 

o 

o 

o 

o 

S-H 

a 

to 

CM 

CO 

u_ 

X 

G 

H 

o 

O 

o 

X 

o 

0 

0 

t-< 

Oh j 

hr 

CO 

Tp 

0 

Q 

< 

LO 

CO 

CM 

G 

T* 

rH 

CO 


o 



G 


Cu 

< 

\ 


\ 

g 

tP 

\ 

cd' 




to 




CO 

t-H 

CM 


00 

o 

00 


CO 

<— i 




jG 

H-H 

Tf 


o 


r-H 

'r 

N 

1 — < 

13 

LO 

j*3 




*G 



‘g 



X 

13 

LO 

h=3 

to 

CM 

tP 

■g 1 

00 

o 

00 


LO 

CN 

CO 



C/5 











C/5 



C/5 

S 




G 




"CJ 




e2 

G 

Ph 

s 

o 

o 

O 


° 

©^ 

o 

E6 

G 

CM 

O 

CM 

1 

G 

d 

o 

O 

o 


o 

o 

o 

£ 

o 


o 



O 

1— J 

LO 

CO 

CO 

CM 

CM 

s 

CO 

o 

CM 

p 

CM 

cn 

u 

! CO 
CO 

CO 

O 

CM 

CM 

rH 


dG 

o 

2 

Tt* 

to 

Tf 

CM 

o 

TP 

s 

CO 

o 


CM 

05 

O 

CJ 


tP 

CM 

oT 


1 





C/5 



O 




CQ 







C/5 

CM 



a 










s 

3 






















G 

Ph 


o 

00 


O 

P 

o 

S 

o 

p 

o 

d 

2 


G 

Q_i 





2 




S 





JG 




O 


d 

d 

d 


d 

d 

d 





p 

p 

CM 


o 

O 

O 


p 

p 

p 


*— J 


1 

1 

| 











r^H 


cd 

cd 

cd 


d 

d 

d 


d 

d 

d 


d 















c 











































G 

Ph 

o 


CO 

p 

p 


p 

p 

p 


p 

p 

r-H 



G 

Oh 


p 


05 


rv _ 

p 

p 


p 

to 


G 

t* 

d 


CO 

i 

LO 

7 


cd 

1 

cd 

1 

7 


1—H 

d 

1 

cd 


G 

JG 

p 

d 


LO 

rd 

d 


d 

CM 

CM 

1 



CM 

td 


G 


05 

p 

oo 


t^ 

p 

p 


to 

Tt 

o 


« 

G 


oq 

P 

p 


r- 







oT 

Ph 



cd 

d 


CM 

d 



CM 

t- 

d 


o 

Ph 



oq 

p 


p 

p 

oo 

CJ 



1 

1 

1 






| 



G 

jG 


id 

d 

d 


cd 

d 

CM 


cd 

id 

CM 

s 















G 

G 




1 










< 3 ; 

to 

G 

Ph 

JG 


p 

p 

CO 


oc 

p 

p 


p 

p 



H 

tg 

G 













*3 


rH 

cd 

d 


CM 

CM 

CO 


cd 

to 

CM 


Q-i 


C5 


p 


p 

p 

CO 


p 

p 

p 

g 

g 

CM 


1 

1 







1 

I 

1 


' G 

CM 


td 

d 

td 


cd 

d 

CO 


TP 

td 

CO 


G 

Ph 

jG 


h- 


_ 


ac 

CO 

LO 


p 

o 

LO 


G 

a 

G 















r-i 

d 

d 


d 

d 

CM 


cd 

o 3 

CM 


<3 

Sg 


O 

p 

OC 1 


t- 

r- 

GO 


05 

p 

p 


LO 


1 

1 

! 



I 



1 

1 

i 



lO 


LO 

td 

TP 


d 

d 

d 


d 

to 

d 


G 

Ph 


CM 

CM 

© 


h- 

LO 

00 


p 

p 

p 



G 

Oh 


P 

00 

LO 


p 

p 



to 

TP 



jG 


CM 

d 

d 


d 

d 

d 


CM 

LO 






CM 

~H 

CM 


CM 

CM 

CM 


d 

o 



o 


t 

1 

t 






1 

1 

1 



o 











I 





jd 

H-J 

G 

Feb. 

|> 

IS 

1 


-G 

HH> 

G 

Feb. 

s 


jG 

G 

-Feb. 

& 

s 





X3 

'r 

-Feb. 

s> 

s 


GG 

'g 

-Feb. 

1 

s 


rG 

Feb. 

H 

k3 




O 

G 

h- 

Nov.- 

1 

d 

£ 


C 

6 

T- 

> 

o 

Z 

A> 

&H 


0 

a 

r- 

Nov.- 

d 

d 





0 

g 

Nov.- 

1 

_d 

ttc 


0 

t^ 

Nov.- 

1 

d 

Uh 


O 

t- 

> 

o 

Z, 

- *q^ 1 


134 




LIMS MAT minus rocket sonde 






Table B15. Statistics for Fort Churchill (59° N, 266° E) 
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